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Gamma-ray bursts (GRBs) have been regarded as standard candles at very high redshift for 
cosmology research. We have proposed a new method to calibrate GRB distance indicators with 
Type la supernova (SNe la) data in a completely cosmology-independent way to avoid the circularity 
problem that had limited the direct use of GRBs to probe cosmology [N. Liang, W. K. Xiao, Y. Liu, 
and S. N. Zhang, Astrophys. J. 685, 354 (2008).]. In this paper, a simple method is provided to 
combine GRB data into the joint observational data analysis to constrain cosmological models; in 
this method those SNe la data points used for calibrating the GRB data are not used to avoid any 
correlation between them. We find that the ACDM model is consistent with the joint data in the 1-er 
confidence region, using the GRB data at high redshift calibrated with the interpolating method, the 
Constitution set of SNe la, the cosmic microwave background radiation from Wilkinson Microwave 
Anisotropy Probe five year observation, the baryonic acoustic oscillation from the spectroscopic 
Sloan Digital Sky Survey Data Release 7 galaxy sample, the x-ray baryon mass fraction in clusters 
of galaxies, and the observational Hubble parameter versus redshift data. Comparing to the joint 
constraints with GRBs and without GRBs, we find that the contribution of GRBs to the joint 
cosmological constraints is a slight shift in the confidence regions of cosmological parameters to 
better enclose the ACDM model. Finally, we reconstruct the acceleration history of the Universe 
up to z > 6 with the distance moduli of SNe la and GRBs and find some features that deviate from 
the ACDM model and seem to favor oscillatory cosmology models; however further investigations 
are needed to better understand the situation. 

PACS numbers: 98.80.Es,98.80.-k,98.80. Jk 



I. INTRODUCTION 



Gamma-ray bursts (GRBs) are the most intense explo- 
sions observed so far and likely to occur in high redshift 
range. Recently, GRB 090423 has been observed at a 
very high redshift, z > 8 The early universe can 
be explored by using GRBs at the high redshift which 
is hardly achievable by Type la supernova (SNe la). In 
recent years, several GRB luminosity relations between 
measurable properties of the prompt gamma-ray emission 
with the luminosity or energy have been proposed as dis- 
tance indicators, many of which are the two-dimensional 
(2D) luminosity relations, such as the isotropic energy 
(I?i so ) - peak spectral energy (Speak) relation (i.e., the 
so-called Amati relation) [2j, the luminosity (L) - spec- 
tral lag ("Hag) relation @, the L - variability (V) relation 
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[H, the L - Speak relation [Bj], the L - minimum rise time 
(tri) relation and the collimation-corrected energy 
(S 7 ) - Speak relation (i.e., the so-called Ghirlanda rela- 
tion) with small scattering [7). In order to reduce the 
large scattering in some of these 2D relations, several 3D 
luminosity relations are also presented [sj-fioj . Liang and 
Zhang first proposed the relation among S; S o, S pea k and 
tb (i.e., the so-called Liang-Zhang relation) [8j, where th 
is the break time of the optical afterglow light curves; 
and Firmani et al. found the relation among L, S pea k 
and To. 45 (i.e., the so-called Firmani relation) [9j, where 
So. 45 is the "high-signal" timescale of the prompt emis- 
sion. Both of the 3D luminosity relations above could 
reduce the scattering in the luminosity or energy versus 
the Speak relation. More recently, Yu et al. find that, for 
the 3D luminosity relations between the luminosity and 
an energy scale S pea k and a timescale (n ag or trt), the 
intrinsic scattering is considerably smaller than that of 
those corresponding 2D luminosity relations [lOj . Other 
GRB luminosity relations have been proposed in many 
lj^ 13 L For reviews of GRB luminosity relations, 



works 
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However, all the luminosity relations of GRBs pre- 
sented above are always empirical, but still without solid 
physical interpretations [16j . Moreover, the luminosity 
relations can not be calibrated with a sufficiently large 
low-redshift GRB sample. Therefore previous calibra- 
tions have been usually obtained by assuming a partic- 
ular cosmological model. Once these GRB luminosity 
relations are calibrated, GRB data could be considered 
as standard candles at very hig h redshift for cosmology 
research [8l fl6U35| (see e.g. [l6l] and [I?} for reviews). In 
2003, Schaefer derived the luminosity distances of nine 
GRBs with known redshifts by using two GRB luminos- 
ity relations {L - jw, L - V) to construct the first GRB 
Hubble diagram [l8f . For 16 GRBs with redshift mea- 
surements, Bloom et al. found a narrow clustering of 
geometrically-corrected Gamma-ray energies within the 
framework of the uniform conical jet model (l9j . Dai 
et al. proposed an approach to consider the Ghirlanda 
relation to constrain cosmological parameters and dark 
energy [20| • Because of the small scatter, the 3D luminos- 
ity relations have been used for cosmological constraint, 
includin g th e Liang- Zhang relation Q, and the Firmani 
relation [24 , |25| . In 2007, Schaefer used five 2D relations 
(L - r lag , L - V, L - E pcak , L - r RT and -E peak ) to con- 
struct the Hubble diagram of 69 GRBs [13]. These GRB 
data have been widely used to constrain cosmology and 
dark energy models in many recent works [28l [29j . and 
joint constraints by combining these GRB data with SNe 
la and the other cosmological probes have been derived 
in [30l - [3"3l | . Instead of a hybrid sample over the whole 
redshift range of GRBs, Takahashi et al. [34J firstly cal- 
ibrated GRB relations (L - T\ ag , L - V) at low redshift 
where distance-redshift relations have been already de- 
termined from SNe la. This method was adopted by 
Bcrtolami and Silva for considering the use of GRBs at 
1.5 < z < 5 calibrated with the bursts at z < 1.5 as 
distance markers to study the generalized Chaplygin gas 
model [35]. All the calibration methods above carried 
out have been derived usually from the ACDM model 
with particular model parameters according to the con- 
cordance cosmology. 

An important point related to the use of GRBs for cos- 
mology is the dependence on the cosmological model in 
the calibration of GRB relations, which had limited the 
direct use of GRBs for cosmology research. In order to 
investigate cosmology, the relations of standard candles 
should be calibrated in a cosmology-independent way; 
otherwise, the circularity problem cannot be avoided eas- 
ily [l6j]. Recently, the possibility of calibrating GRBs in 
a low-dispersion in redshift near a fiducial redshift has 
been proposed [H, H3], which has been developed fur- 
ther based on Bayesian theory (38|. However, the GRB 
sample available now is far from what is needed to cali- 
brate the relations in this way [ill [38j] . Many of the works 
treat the circularity problem with statistical approaches 
[HI, [23|, l39l - [4l| . A simultaneous fit of the parameters in 
the calibration curves and the cosmology is carried out 
to find the optimal GRB relation and the optimal cos- 



mological model in the sense of a minimum scattering in 
both the luminosity relations and the Hubble diagram 
[l8j . Firmani et al. proposed a Bayesian method to get 
around the circularity problem [23j . Li et al. presented a 
global fitting analysis for the Ghirlanda relation to deal 
with the problem [39]. Amati et al. used the _E; so - 
-Epcak relation to measure the cosmological parameter by 
adopting a maximum likelihood approach [401 ] . More re- 
cently, Wang has shown that the current GRB data can 
be summarized by a set of model independent distance 
measurements, with negligi ble loss of information [4l| . 
which is followed by |42H44j . However, an input cosmo- 
logical model is still required in doing the joint fitting, the 
circularity problem can not be circumvented completely 
by means of these statistical approaches. 

In our previous paper (45[ , we presented a new method 
to calibrate several GRB luminosity relations in a com- 
pletely cosmology-independent manner. Our method 
avoids the circularity problem more clearly than previous 
cosmology-dependent calibration methods. It is obvious 
that objects at the same redshift should have the same 
luminosity distance in any cosmology. Therefore, the lu- 
minosity distance at any redshift in the redshift range 
of SNe la can be obtained by interpolating (or by other 
mathematical approach) directly from the SNe la Hubble 
diagram. Using the interpolation method, we calibrated 
seven GRB luminosity relations (L - ri ag , L - V, L - 

Epeak, L - Trt, Ej - -E pca k, -Esq " -Epeak, and Ei so -E pca k 



tb) with the data compiled by Schaefer [17J. Then if fur- 
ther assuming these calibrated GRB relations valid for 
all long GRB data, we can use the standard Hubble di- 
agram method to constrain the cosmological parameters 
from the GRB data at high redshift obtained by utilizing 
the relations [45| . These distance data of GRBs are so 
far the most cosmology-independent GRB distance indi- 
cators. Following this cosmology-independent GRB cali- 
bration method from SN la [45[ , the derived GRB data at 
high redshift range can be used to constrain cosmologi- 
cal models without circularity problem [45Tj47i l49j - [5l| . 
Capozziello and Izzo first used the Eiso-Ep^-t^ relation 
and the E 1 - E pcak relation calibrated with the so-called 
Liang method to derive the related cosmography param- 
eters [deceleration parameters (q), jerk parameters (j), 
and snap parameters (s) [481]] . which are only related to 
the derivatives of the scale factor without any a priori 
assumption [46[ . Wei and Zhang used the Amati relation 
calibrated with the interpolation method to reconstruct 
the acceleration history of the Universe and constrain 
cosmological models [49|, Ho| • 

Furthermore, we proposed another approach to cali- 
brate GRB relations by using an iterative procedure [HJ , 
which is a nonparametric method in a model indepen- 
dent manner to reconstruct the luminosity distance at 
any redshift in the redshift range of SNe la [531-55], Sim- 
ilar to the interpolation method, Cardone et al. con- 
structed an updated GRBs Hubble diagram on six 2D 
relations calibrated by local regression from SNe la [56| . 
This method and GRB data have been used to cosmo- 
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logical constraints with SNe la in some following works 
[ia . Kodama et al. presented that the L - i? pca k re- 
lation can be calibrated with the empirical formula fitted 
from the luminosity distance of SNe la [58| . This method 
has been used to constrain cosmological parameters by 
combining these GRB data with SNe la in a following 
work [59j]. However, it is noted that this calibration pro- 
cedure depends seriously on the choice of the formula and 
various possible formulas can be fitted from the SNe la 
data that could give different calibration results of GRBs. 
As the cosmological constraints from GRBs are sensitive 
to GRBs calibration results [41| . the reliability of this 
method should be tested carefully. Moreover, as pointed 
out in 41], the GRB luminosity relations which are cali- 
brated by this way are no longer completely independent 
of all the SNe la data points; therefore these GRB data 
can not be used to directly combine with the whole SNe 
la dataset to constrain cosmological parameters and dark 
energy. 

It is easy to find that the number of SNe la data points 
that have been used in the linear interpolating procedure 
to obtain the GRB data is relatively small compared to 
the whole SNe la sample. In order to combine GRB data 
obtained by our interpolation method into the joint ob- 
servational data analysis to constrain cosmological mod- 
els, we can exclude the SNe points that have been used in 
the interpolating procedure from the SNe la sample used 
to the joint constraints. In this work, with the updated 
distance moduli of the 42 GRBs at z > 1.4 obtained by 
the interpolating method [45[ from the Union set of 307 
SNe la [6J], we provide a simple method to ensure that 
the calibrated GRB data are independent of the SNe la 
data used in the joint data fitting to constrain cosmolog- 
ical models, by using the Constitution set of 397 SNe la 
65] , along with the cosmic microwave background radia- 
tion (CMB) observation from the five-year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP5) result 1661 . (R, l a , 
ri{,/i 2 ); the baryonic acoustic oscillation (BAO) [67] obser- 
vation from the spectroscopic Sloan Digital Sky Survey 
(SDSS) Data Release 7 (DR7) galaxy sample (d . 2 , 
do.35)', the 26 baryon mass fraction in clusters of galaxies 
(CBF) from the x-ray gas observation (/ gas ) HH; and the 
11 Hubble parameter versus redshift data, H(z) [z3,[z3. 
Our goal is to determine the contribution of GRBs to the 
joint cosmological constraints in the confidence regions of 
cosmological parameters by comparing to the joint con- 
straints with GRBs and without GRBs. Finally we also 
reconstruct the acceleration history of the Universe up to 
z > 6 with the distance moduli of SNe la and GRBs. 



II. OBSERVATIONAL DATA ANALYSIS 

The updated distance moduli of the 42 GRBs at z > 
1.4 are obtained by the five GRB luminosity relations 
that are calibrated from the 27 GRBs at z < 1.4 with the 
interpolating method [45[ using the Union set of 307 SNe 
la [64l ] . The Union set includes the Supernova Legacy 



Survey (SNLS) ^ and the ESSENCE Survey J6l|, the 
extended dataset of distant SNe la observed with the 
Hubble Space Telescope [62], and the formerly observed 
SNe la data. For the SNe data in joint analysis, we use 
the new Constitution set of 397 SNe la [65[ , which com- 
bine 90 CfA3 SNe la sample with the Union set [Hj]. We 
also consider the shift parameters set (R, l a , D,i,h 2 ) [66| . 
the BAO parameters set (d . 2 , do.35) from SDSS DR7 
[68j . the 26 / gas data from the x-ray gas observation in 
clusters of galaxies 69], and the 11 H{z) data from the 
differential ages of passively evolving galaxies [7(3] and 
from the BAO peak position [7l| . 

In our previous paper [45L we used the data for 69 
GRBs compiled by Schaefer pjj and adopted the data 
for 192 SNe la compiled by Davis et al. [63| to calibrate 
seven GRB luminosity relations {L - ri ag , L - V, L - 

-Speak, L - Trt, Ej - -Ep Ca k, E[ so - Speak) and -Eiso-£pc a k- 

tb) at z < 1.4 by using the interpolation method. In this 
work, we use the Union set of 307 SNe la [64| to deter- 
mine the values of the intercepts and the slopes of GRBs 
luminosity relations calibrated with the GRB sample at 
z < 1.4 by using the linear interpolation method. We do 
not include the 90 SNe la data from CfA3 [65] due to 
their extremely low redshift (z < 0.1), which would not 
affect the calibrated results for GRB luminosity relations 
at 0.17 < z < 1.4. The 2D luminosity relation of GRBs 
can be generally written in the form 

logy = a + blogx, (1) 

where y is the luminosity (L/erg s -1 ) or energy (E~ ( / 
erg); x is the GRB parameters measured in the rest 
frame, e.g., r lag (l + z)-V(0.1 s), U(l + z)/0.02, £ poak (l + 
z)/(300 keV), trt(1 + z) _1 /(0.1 s), for the corresponding 
2D relations. We adopt the data for these quantities from 
Ref. [l7| including five 2D GRB luminosity relations (L 
- riag, L - V, L - -Epcak, L - r RT and E 1 - -E pcak ). For 
the E~ - -E pea k relation, in order to calculate the total 
collimation-corrected energy Ej, one needs to know the 
beaming factor, -Fbc a m = (1 — cos6>j G t), where the value 
of the jet opening angle 0j Ot is related to the jet break 
time (tb) and the isotropic energy for an Earth-facing jet, 
= £ ; 7,iso/10 52 erg [72j]. When calculating E^ t - lso , 
we also use the interpolation method from SNe la to avoid 
the circularity problem. We determine the values of the 
intercept (a) and the slope (b) with their 1-cr uncertain- 
ties calibrated with the GRB sample at z < 1.4 by using 
the linear interpolation methods from the Union set of 
307 SNe la. 

Further assuming that these GRB luminosity relations 
do not evolve with redshift, we are able to obtain the 
luminosity (L) or energy (-Ely) of each burst at high 
redshift (z > 1.4). Therefore, the luminosity distance 
{d£) can be derived. A distance modulus can be cal- 
culated as fx = 51og(c?i/Mpc) + 25. The weighted av- 
erage distance modulus from the five relations for each 
GRB is fi = (J2if i i/ cr fi i )/(J2i' T ^i 2 )^ wrt h its uncertainty 
a ^ = tT ^ i 2 ) _1//2 i where the summations run from 1 to 
5 over the five relations with available data. For more de- 
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FIG. 1: Hubble Diagram of 397 SNe la (red dots) and the 
69 GRBs (circles) obtained using the interpolation methods. 
The 27 GRBs at z < 1.4 are obtained by linear interpolating 
from SNe la data (black circles), and the 42 GRBs at z > 1.4 
(blue circles) are obtained with the five relations calibrated 
with the sample at z < 1.4 using the interpolation method. 
The curve is the theoretical distance modulus in the concor- 
dance model (w = — 1, Quo = 0.27, JIa = 0.73), and the 
vertical dotted line represents z = 1.4. 



tails of the calculation, see [17j and 1451 . We have plotted 
the Hubble diagram of 397 SNe la [6fJ which combine 90 
CfA3 SNe la sample with the Union set [gij and the 69 
GRBs obtained using the interpolation methods in Fig 1. 
The distance moduli of the 27 GRBs at z < 1.4 are ob- 
tained by using the linear interpolation method directly 
from the Union SNe data. The 42 GRB data at z > 1.4 
are obtained by utilizing the five relations calibrated with 
the sample at z < 1.4 using the interpolation method. 

Constraints from the SNe la data and the GRB data 
can be obtained by fitting the distance moduli fJ.(z). A 
distance modulus can be calculated as 



M = 5 log 



Mpc 



25 = 51og 10 Dl — Mo, 



(2) 



where Hq = 5 log 10 h + 42.38, h = i? /(100km/s/Mpc), 
Hq is the Hubble constant, and the luminosity distance 
D L can be calculated by 



D L = H d L = (1 + z)fl k 1/2 sinn 



1/2 



dz' 



E(z>) 



(3) 



where fl^ = 1 — SImo — ^de, sinn(a;) is sinh for J7k > 0, 
sin for D,^ < 0, and x for = 0; and E(z) can be given 
by 



-no 



(4) 



= [0 M o(l + -z) 3 + f^DEexp 
+fi k (l + z) 2 ] 1 / 2 , 



:W Z1 + W( f ) ^ 
1 + z' 



where the equation of state (EoS) of dark energy w(z) 
is determined by the choice of the specific dark energy 
model. The % 2 value of the observed distance moduli 
can be calculated by 



N 

£ 

i=l 



(5) 



where /i bs(^i) is the observed distance modulus for 
the SNe la and/or GRBs at redshift Zi with its er- 
ror (Juj is the theoretical value of distance mod- 
ulus from a dark energy model which can be calcu- 
lated from Eqs. -(4). The nuisance parameter /xo 
(or h) can be obtained by marginalizing the likeli- 
hood functions of over /iq (or h) for all values [2(| 
or treated by following an effective approach [73| to 
expand the x 2 to xtif^o) = ^Mo — 2£>,uo + C, with 
^4 = EV<, B = E^obsN-Slog^^]/^, and 
C = E [Mobs(zi) - 5 log 10 Di] 2 /cr 2 .. Thus the x 2 has a 
minimum value for /i = £>/A at x 2 which given by 



B 2 



(6) 



Therefore we can minimize x 2 instead of minimizing x 2 . 

For the CMB observation from the WMAP5 results 
[66| . the two shift parameters R and i , together with 
the baryon density fraction of the Universe (cj& = Vl^ti 2 ) 
can provide an efficient summary of CMB data to con- 
strain cosmological models. The shift parameter R can 
be expressed as 



R = £\io^k smn 



1/2 



(7) 



o 

is the redshift of recombination which is given 



where z rec 

by [z3l, Zrcc = 1048[l+0.00124w 



-0.738 



(i+ 5l (Wi 2 ) 92 )], 



91 



0.0783a;r; - 238 (l-|-39.5w ; 7 a763 )- 1 and ff2 = 0.560(1- 



21.1wJ- 81 )- 



The shift parameter ^ Q can be expressed as 



J a = 7T 



^^sinn^ 2 

(-^rcc) 



(8) 



where r s (z rcc ) is the comoving sound horizon at photo- 
decoupling epoch, 



f s (-^rec) 



1 

So 

Cl lcc 



c s (z) 

E( z y 



(9) 



1/2 
M0 



fL/l 2 



1-1/2 



da 



with the sound speed c s = (l/v / 3)[l 
a(31500wfe(T'cMB/2.7K)- 4 )]- 1 / 2 0, and T CM b 
2.725 K [66j]. From the WMAP5 measurement,!] 
best-fit values of (R, l a , 100a>b) for a flat prior are |66l ] 



CMB 




1.710 ±0.019 
302.10 ±0.86 
2.2765 ±0.0596 



(10) 



5 



The x 2 value of the CMB observation from WMAP5 can 
be expressed as 



XcMB 



— AP^a/trC 



CMB^CMB 



X AP 



CMB, 



where 



AP 



CMB 



R- 1.710 
l a - 302.10 
100w b - 2.2765 



(11) 



(12) 



and the corresponding inverse covariance matrix is [661 ] 

2809.73 -0.133381 158.356 \ 
• r\ 1L; - j -0.133381 2.21908 19.7195 . (13) 
158.356 19.7195 465.728 / 

It is noted that we use the priors following [76j, fib = 
0.022765//i 2 and h = 0.705 [66j], when calculating the 
value of Xcmb- 

For the BAO observation [63, from the SDSS DR7 
galaxy sample [68[ , we use the measurement of d z at z = 
0.2 and z = 0.35, where d z can be expressed as (68[ 



d z = 



r s (zd) 



-Dv(zbao) 



(14) 



where z& is the drag epoch at which baryons were 
released from photons which is given by [751 ] . Zd = 

1291(ft M o^)°- 251 [(l + 6i^ b2 )]/[l + 0.659(^Mo^ 2 ) a828 ], 
h = 0.313(fi M ofr 2 r°- 419 [l + 0.607(n M0 h 2 )°- 67i }-\ and 



6 2 = 0.238(fW» 2 ) ' 223 , and D v can be given by [67] 

f ZBAO dz 
'o 



Dv{zbao) 



1 



zbao 



H l E(zbao) 



\ 2 1 1/3 



E{ z y 



(15) 



From the SDSS DR7 measurement, the best- fit values are 



BAO = 



do. 2 
^0. 35 



0.1905 ±0.0061 
0.1097 ±0.0036 



(16) 



The x 2 value of the BAO observation from SDSS DR7 
can be expressed as [o8| 



Xbao 



APbaqCbao 1 AP B ao, 



(17) 



where 



AP 



BAO 



d .2 - 0.1905 
do.35- 0.1097 



(18) 



and the corresponding inverse covariance matrix is [68 



'BAO 



30124 -17227 \ 
-17227 86977 ) 



(19) 



The baryon mass fraction in clusters of galaxies from 
the x-ray gas observation can be used to constrain cosmo- 
logical parameters on the assumption that the gas mass 
fraction in clusters is a constant and thus independent 



of redshift. The baryon gas mass fraction / gas can be 
presented as [69| 



jSCDM 



f ( z ) - \\^A Wl2/3 

/gasl j_A[ d A (z) 1 ' 



(20) 



where A = [6fi b (2/i) 3/2 ]/[(l + a)fi M o], a = 0.19Vh, b is a 
bias factor motivated by gas dynamical simulations, and 
dA = di/(l + z) 2 is the theoretical value of the angular 
diameter distance from cosmological models, d^ CDM is 
the angular diameter distance corresponding to the stan- 
dard cold dark matter (SCDM) universe (f^Mo = 1 for a 
flat universe). Here we adopt the usually used 26 obser- 
vational /ga S data [g9j] to constrain cosmological models. 
The x 2 value of cluster's baryon gas mass fraction is 



N=26 rj-obs 



XCBF 



E 

i=l 



[/gas ( Z i) /gas( 2 



(21) 



Following [77J, |78| , we treat A as a nuisance parameter to 
expand the Xcbf to Xcbf(^) = ^4A 2 -2BA + C, with A — 
E[/gas,i/o-/ gas ,i] 2 , B = E[/gas,i/gas,i]/a-^ asii , and C = 
E [W<^as,d 2 , where / gaSil = ]4 CDM (z)/d A (z)] 2 / 3 . 
Thus the Xcbf a minimum value at Xcbf wmcn is 
given by 



y 2 -C B2 
Xcbf — ° ^ 



(22) 



The Hubble parameter H(z) can be derived from the 
derivative of redshift with respect to the cosmic time, 



H(z) = 



1 dz 
1 + z dt ' 



(23) 



From the Gemini Deep Deep Survey (GDDS) [79] obser- 
vations of differential ages of passively evolving galaxies 
and other archival data [80(, the H(z) data at nine differ- 
ent redshifts (0.09 < z < 1.75) have been obtained [70| . 
Recently, using the BAO peak position as a standard 
ruler in the radial direction, H = 83.2±2.1km/s/Mpc at 
z = 0.24, and H = 90.3 ± 2.5km/s/Mpc at z = 0.43 have 
been obtained [7l|. To constrain cosmological models, 
the x 2 value of the 11 H(z) data is 



N S A 1 [H ohs (zi) - H(zi)] 



i=l 



(24) 



'H,; 



The nuisance parameter Hq is also marginalized follo wing 
the procedure used in calculating x 2 and Xcbf 



III. CONSTRAINTS ON COSMOLOGICAL 
MODELS FROM THE JOINT DATA WITH GRBS 

The distance moduli of the 42 GRBs at z > 1.4 are 
obtained by the five GRB luminosity relations that are 
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calibrated from the 27 GRBs at z < 1.4 with the in- 
terpolating method using directly the distance moduli of 
adjacent SNe la. In the interpolating procedure to obtain 
the distance moduli of the 27 GRBs at z < 1.4, we have 
used only 40 SNe la data points from the Union set of 
307 SNe la. In order to combine GRB data into the joint 
observational data analysis to constrain the dark energy 
models, we exclude the 40 SNe points from the SNe la 
sample used to the joint constrains. Therefore the re- 
maining SNe la data points are completely independent 
of the distance moduli of the 42 GRBs at z > 1.4. Those 
excluded 40 SNe are listed in the Appendix. Since the 
reduced 357 SNe la, 42 GRBs, CMB, BAO, as well as 
CBF and H(z) are all effectively independent, we can 
combine the results by simply multiplying the likelihood 
functions. Thus the cosmological parameters can be fit- 
ted with the combined observable data by the minimum 
X 2 method. The total x 2 with the SNe + GRBs + CMB 
+ BAO + CBF + H{z) dataset is 

X 2 = X^,{SN+GRB} + XCMB + XbAO + XcBF + Xh ■ (25) 

The best-fit values for these parameters can be deter- 
mined by minimizing the total % 2 . For comparison, SNe 
+ CMB + BAO + CBF + H(z) without GRBs have 
been used to show the contribution of GRBs to the joint 
cosmological constraints, and we also consider the joint 
constraints with 397 SNe + CMB + BAO + CBF + 
H(z). In addition, some different data set such as SNe + 
GRBs + CMB + BAO, SNe + GRBs + CMB, and SNe 
+ GRBs + BAO have also been used in the cosmological 
constraints. 

We use the Akaike information criterion (AIC) [8l| and 
the Bayesian information criterion (BIC, the so-called 
Schwarz information criterion) [82[ to select the best-fit 
models. Liddle examined the use of information crite- 
ria in the context of cosmological observations [83[ . The 
AIC is defined as HH AIC = -21nJz? max + 2k, where 
-S^max is the maximum likelihood, and k the number of 
parameters of the model. Models with too few param- 
eters give a poor fit to the data and hence have a low 
log-likelihood, while those with too many are penalized 
by the second term. The best model is the model which 
minimizes the AIC 83] . The BIC can be defined as [82[ 
BIC = — 21nJ§f max + klnN, where N is the number of 
data points used in the fit. Note that for Gaussian errors, 
Xinin = — 2 In .5f m ax, and the difference in BIC can be sim- 
plified to BIC = Ax^in + Afcln N. The AIC gives results 
similar to the BIC approach, although the AIC is not 
strict enough on models with extra parameters for any 
reasonably sized data set (In A > 2). Therefore, for com- 
paring cosmological models from the joint data set, we 
only compare ABIC measured with respect to the best 
model. A difference in BIC of 2 is considered positive 
evidence against the model with the higher BIC, while a 
ABIC of 6 is considered strong evidence [HI . 

The combined data sets are used to constrain cosmo- 
logical parameters and dark energy. Here we consider 
three cosmological models, the ACDM model with dark 



energy EoS w = — 1, the u>CDM model with a constant 
EoS, and the Chevallier-Polarski-Lindcr (CPL) model in 
which dark energy with a parameterization EoS as [84| 



w(z) = w + w a (l - a) = w a + w a 
For the ACDM model, Eq. (J3J becomes 



1 + 2 



1/2 



(26) 



(27) 



where J7mo + + f2k = 1< For the wCDM model with 
a constant EoS for a flat universe prior, 

E(z) = [fi M0 (l + zf + (1 - QmoXI + z) 3(1+Wo) ] (28) 

and considering w(z) as CPL parameterization model for 
a flat universe prior, Eq. (|4]) becomes 



E(z) 



Omo(1 + zf + (1 - ^mo) 

(1 + z )3(l+wo+Wa) e -3w aJ ^; 



(29) 



1/2 



In Fig. 2, we show the joint confidence regions in the 
(^ M0 ,^a) plane for the ACDM model. With 357SNe 
+ GRBs + CMB + BAO + CBF + H(z), the 1-cr 
confidence region for (J7mo,^a) of the ACDM model 
are (Q M oM = (0.2751°;°^, 0.723±°;°i£), with x 2 min = 
494.476 for 438 degrees of freedom. For comparison, 
fitting results from the joint data without GRBs are 
also given in Fig. 2. With 357SNe + CMB + BAO 
+ CBF + H(z), the best-fit values are (Omoj^a) = 
(Q.270±g;g^, 0.730±g;|££) with xL in = 449.001 for 396 de- 
grees of freedom. We present the best-fit value of Omo, 
f^A with lcr uncertainties, and Xmin' Xmin/^ ^ as weu as 
AIC, BIC for the ACDM model in Table 1. 

From comparing to the joint constraints with GRBs 
and without GRBs, we can see that the contribution of 
GRBs to the joint cosmological constraints is a slight 
shift between the best-fit values near the line that rep- 
resents a flat universe, towards a higher matter density 
Universe by AJImo = 0.005, compared to the joint con- 
straints without GRBs. It is noted that the obtained 
errors for these parameters with and without GRBs are 
essentially unchanged, because the number of SNe la data 
points dropped is similar to the number of GRBs when 
including GRBs in the joint fitting. This model has the 
lowest BIC compared to other models tested from the 
joint data [357SNe + CMB + BAO + CBF + H(z)} 
with GRBs, so ABIC are measured with respect to this 
model. 



Figure 3 shows the joint confidence regions in the 
(fiM0;w) plane for the wCDM model with a constant 
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ACDM model (with GRBs) 


ACDM model (without GRBs) 


wCBM model (with GRBs) 


wCDM model (without GRBs) 


Omo 

n A 

w 


n 79 o + 0.017 

tu = — 1 
494.476 


270+ 016 
u. rdu_ 017 
w = — 1 
449.001 


269+ 013 

U.ZD»_ Q Q14 

Oa = 1 — Hmo 
-0 99+ 07 

u - J -0.07 

494.532 


267+ 013 

Qa = 1 — ^MO 

-0 98+ 07 
448.808 


xL„/ d ° f 


1.129 


1.134 


1.129 


1.133 


AIC 


500.476 


454.001 


500.532 


454.808 


BIC 


512.736 


466.960 


512.792 


466.767 



TABLE I: The best-fit value of Qmo, ^a, and w with la uncertainties, and Xmim Xmm/dof, as well as AIC, BIC 
for the ACDM model and for the flat tuCDM model with SNe+GRBs+CMB+BAO+CBF+#(z) ( w ; t h GRBs) and 
SNc+CMB+BAO+CBF+flYz) (without GRBs). 
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FIG. 2: The joint confidence regions in the 
(Omo^a) plane for the ACDM model with 
SNe+GRBs+CMB+BAO+CBF+#(,z). The contours corre- 
spond to l-<7 and 2-cr confidence regions, and the red point 
is the best-fit value [SNe+GRBs+CMB+BAO+CBF+tf^)]. 
The results for the cases with and without GRBs are indi- 
cated by the solid lines and the dot-dashed lines, respectively. 
The blue line represents a flat universe (JImo + Oa = 1). 
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FIG. 3: The joint confidence regions in the (Amo,w) plane 
for the dark energy model with a constant w for a flat 
universe with SNe+GRBs+CMB+BAO+CBF+#(z). 
The contours correspond to 1-er and 2-a confi- 
dence regions, and the red point is the best-fit value 
[SNc+GRBs+CMB+BAO+CBF+H(z)]. The results for the 
cases with and without GRBs are indicated by the solid lines 
and the dot-dashed lines, respectively. 



EoS for a flat universe. With 357SNe + GRBs + 
CMB + BAO + CBF + H(z), the 1-cr confidence region 
for (tt-M_o,w) of the flat roCDM model are (flM,w) = 
(0.269i°;°i|, -0.99tg$), with Xmin = 494.532 for 438 
degrees of freedom. For comparison, fitting results from 
the joint data without GRBs are also given in Fig. 3. 
With 357SNe + CMB + BAO + CBF + H{z), the best- 
fit values are {n M0 ,w) = (0.267^0^^,-0.98^;^), with 
Xmin — 448.808 for 396 degrees of freedom. We present 
the best-fit value of JImoj with la uncertainties, and 

Xmin: Xmin/ dof > aS wel1 aS AIC > BIC for the nat »CDM 

model in Table 1. 

Comparing to the joint constraints with GRBs and 



without GRBs, we can see that the contribution of GRBs 
to the joint cosmological constraints is a slight shift that 
adding the best-fit value of fijyio to 0.002, and subtract- 
ing the best-fit value of w to —0.01 to enclose the ACDM 
model (w = -1). From the joint data (357SNc + CMB 
+ BAO + CBF + H(z)) with GRBs for the flat wCDM 
model, we obtain ABIC = 0.056 with respect to the 
ACDM model. This has the lowest BIC compared to 
other models tested from the joint data [357SNe + CMB 
+ BAO + CBF + H{z)] without GRBs, so ABIC are 
measured with respect to this model. Therefore, for the 
ACDM from the joint data without GRBs , we obtain 
ABIC = 0.193 with respect to the flat wCDM model. 
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Qmo 


0.269 


0.269 


0.269 


0.269 


0.263 


0.307 




494.530 


448.726 


504.496 


458.857 


455.554 


452.839 




1.132 


1.136 


1.160 


1.147 


1.145 


1.138 


AIC 


502.530 


456.726 


512.496 


464.857 


463.554 


460.839 


BIC 


518.877 


472.672 


528.825 


482.853 


479.530 


476.815 



TABLE II: The best-fit value of the corresponding (wo, w a ) with la uncertainties and the best-fit value of Qmo, and Xmim 
xLn/dof, as well as AIC, BIC for the CPL model in a flat universe, with 357SNe+GRBs+others(CMB+BAO+CBF-|-fl'), 
357SNe+others (without GRBs), 397SNe+others, and with 357SNe+GRBs+C(CMB)+B(BAO), 357SNe+GRBs+CMB, and 
357SNe+GRBs+BAO. 
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Wo 

FIG. 4: The joint confidence regions in the 
(wo,w a ) plane for the CPL model in a flat universe 
with 357SNe+GRBs+others(CMB+BAO+CBF+#) , 

357SNe+others (without GRBs), and 397SNe+others. 
The contours correspond to 1-a and 2-a con- 
fidence regions, and the red point is the best- 
fit value (357SNe+GRBs+CMB+BAO+CBF+#). 
The solid lines represent the results of 
357SNe+GRBs+others(CMB+BAO+CBF+i7) . The 
dash-dotted lines represent the results of 357SNe+others 
(without GRBs). The dashed lines represent the results of 
397SNe+others. 



For the flat CPL model, we find that the best-fit pa- 
rameters with 357SNe + GRBs + CMB + BAO+ CBF 
+ H(z) are (w ,w a ,n M o) = (-0.98,-0.02,0.269), with 
Xmin — 494.53 for 437 degrees of freedom. Figure 4 
shows the joint confidence regions in the (wo,w a ) plane, 
while fixing J7mo = 0.269. For comparison, fitting results 
from the joint data without GRBs are also given in Fig. 
4. With 357SNe + CMB + BAO + CBF + H(z), the 




-1.2 -1.0 -0.8 -0.6 



wo 

FIG. 5: The joint confidence regions in the (wo,w a ) 
plane for the CPL model in a flat universe with 
357SNe+GRBs+CMB+BAO, 357SNe+GRBs+CMB and 
357SNe+GRBs+BAO. The contours correspond to 1-a and 
2-a confidence regions, and the red point is the best-fit value 
(357SNe+GRBs+CMB+BAO). The solid lines represent the 
results of 357SNe+GRBs+CMB+BAO. The dashed lines 
represent the results of 357SNe+GRBs+CMB. The dash- 
dotted lines represent the results of 357SNe+GRBs+BAO. 
The straight line near the top is given by w a + wo = 0. 



best-fit values are (w ,w a , Qmo) = (—0.95,-0.14,0.269), 
with x^in = 448.726/395. Fitting results with 397SNe 
+ CMB + BAO + CBF + H(z) are also given in 
Fig. 4 and the best- fit values are (wq, w a , ^mo) = 
(-0.99,0.00,0.269), with x min = 504.496/435. In ad- 
dition, fitting results from the joint data with 357SNe + 
GRBs + CMB + BAO, 357SNe + GRBs + BAO, and 
357SNe + GRBs + CMB are given in Fig 5. We note that 
the contours for 357SNe + GRBs + BAO are sharply 
cutoff at the top, near the line given by wq + w a = 0, 
as shown in the figure. This is due to w(z > 1) < or 
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u'o + w a < in the CPL parameterization (i.e, the early 
universe is matter dominated) implicitly required by the 
BAO data; this was also noted by Kowalski et al 
We present the best-fit value of the derived wq, w a with 
la uncertainties, and the best-fit value of SImo, Xmin> 
X^ in /dof, as well as AIC, BIC with differently combined 
data set in Table 2. 

Comparing to the joint constraints with GRBs and 
without GRBs, we can see that the contribution of 
GRBs to the joint cosmological constraints is a slight 
shift in the (wo,w a ) plane to enclose the ACDM model 
(wq = — l,w a = — 0). From the joint data [357SNe 
+ CMB + BAO + CBF + H{z)] with GRBs for the 
fiat CPL model, we obtain ABIC = 6.141 with respect 
to the ACDM model, indicating a strong preference for 
the ACDM model; while from the joint data without 
GRBs for the flat CPL model, we obtain ABIC = 5.905 
with respect to the flat wCDM model. Meanwhile com- 
paring to Fig. 4 and 5, we find the effects of adding 
the 26 /g as data and the 11 H(z) data are not very 
significant, suggesting that the model parameters are 
strongly constrained by SNe+GRBs+CMB+BAO. Com- 
paring to the joint constraints with SNe+GRBs+CMB 
and SNe+GRBs+BAO, we find the confidence regions 
with SNe+GRBs+CMB seem to close to the joint con- 
straints with SNe+GRBs+CMB+BAO. It indicates that 
the contribution of CMB data to the joint cosmological 
constraints (SNe+GRBs) is more significant compared to 
that of BAO data. 

From Fig. 2-5 and Table 1 - 2, we can find that the 
ACDM model is consistent with the joint data in the 1- 
<7 confidence region. Comparing to the joint constraints 
with GRBs and without GRBs, we can find the effect 
of GRBs to the joint cosmological constraints, although 
the contribution of GRBs to the cosmological constraints 
would not be sufficiently significant compared to that of 
SNe la at present. This is mainly caused by the relatively 
large statistical scattering in the GRB relations and the 
relatively small data set of GRBs compared to that of 
SNe la currently. 



IV. RECONSTRUCTING THE ACCELERATION 
HISTORY OF THE UNIVERSE 



where f(z) = Di,{z)/h, f{z) s n represents the smoothed 
luminosity distance at any redshift z after the n th it- 
eration, f(z)g denotes a guess background model, and 
f obs (zi) is the observed one from the SNe la data, as 
well as the normalization parameter N(z)~ 1 = Ni(z), 
Ni{z) =exp[- (ln 2 ((l + z)/(l + z l )))/(2A 2 )]. It has 
been shown that the results are not sensitive to the cho- 
sen value of A and the assumed initial guess model. Here 
we use a wCDM model with w — —0.9 and Omo = 0.28 
as the guessed background model and we choose A = 0.6 

The best iterative result is obtained by minimizing 

Xl = j2(^i)n - »° bS (Zi)) 2 rt obs y (31) 
i 

Once the \n reaches its minimum value, we stop the it- 
erative process and get the best result (Xmin)' w ^h the 
1-cr uncertainties corresponding to x 2 = Xmin + 1 ■ Figure 
6 shows the computed Xn f° r the reconstructed results 
at each iteration for the SN la and GRB data. We find 
that when at n = 16 a minimum value of Xn is obtained. 




20 40 60 80 100 



Iteration number 

FIG. 6: Computed Xn f° r the reconstructed results at each 
iteration for the SN la and GRB data. The minimum value of 
Xn is obtained at n — 16, and the l-a uncertainties of recon- 
structed results can be obtained at n — 9,41, corresponding 
to Ax 2 = 1, respectively. 



Following a well-known procedure in the analysis of 
large scale structure, Shafieloo et al used a Gaussian 
smoothing function rather than the top hat smoothing 
function to smooth the noise of the SNe la data directly 
[53L l54j . In order to obtain important information on in- 
teresting cosmological parameters expediently, when do- 
ing the Gaussian smoothing, lndi(z), rather than the 
luminosity distance di,(z) or distance modulus fJ,(z), is 
studied by an iterative method [53l . [HJ . Here we follow 
the iterative procedure and adopt results from Ref. [55[ , 



By using the iterative approach to reconstruct the lu- 
minosity distance at any redshift in the redshift range 
of SNe la, we have calibrated GRB luminosity relations 
in a model independent manner [52| . Here we recon- 
struct the Hubble parameter H(z), the deceleration pa- 
rameter q(z), and the EoS of dark energy w{z) from the 
best iterative result of f(z) with the distance moduli of 
SNe la and GRBs obtained by the interpolating method. 
The Hubble parameter can be give n by differentiating the 
smoothed luminosity distance [55| . 



In /(*)' = In + N(z) V [ In f° bs ( Zi ) ~ In /(z^-i] N^O) 



Hi*) = -jz 



d_ 

dz 



(1 + *) 



(32) 
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which contains the information on Ho. Then the deceler- 
ation parameter q(z) of the expanding universe and the 
EoS of dark energy can be obtained |55| , 



q(z) = (1+Z) 



H'(z) 
H(z) 



I: 



w{z) — 



1 



-i + |(i + z)-gyg 



(33) 



(34) 



where the prime (/) denotes the derivative with respect 
to z. 



N 

X 




FIG. 7: The reconstructed H(z) with the 1-a uncertainties. 
The red line indicates the best-fit value of the reconstructed 
H(z) which is obtained at n — 16 and the blue lines represent 
the fit values of the reconstructed H(z) within 1-a uncertain- 
ties at n = 9, 14, 19, 24, 29, 34, 39, 41. The dashed line is the 
theoretical values of the ACDM model. 



The reconstructed H (z) with the 1-a uncertainties are 
shown in Fig. 7. For comparison, the theoretical values 
of the ACDM model are also shown in Fig. 7. We can 
find that there is a deviation between the reconstructed 
H(z) and the theoretical values of the ACDM model with 
the 1-a uncertainty at 1 < z < 4, where the GRB data 
clearly dominate. The reconstructed q{z) with the 1-a 
uncertainties are shown in Fig. 8. We can find that the 
transition redshift at which the expansion of the universe 
from deceleration (q(z) > 0) to acceleration (q(z) < 0) 
is zt — 0.38^0 03: & t relatively low redshifts, where the 




FIG. 8: The reconstructed q(z) with the 1-a uncertainties. 
The red line indicates the best-fit value of the reconstructed 
q(z) which is obtained at n = 16 and the blue lines represent 
the fit values of the reconstructed q(z) within 1-a uncertain- 
ties at n = 9, 14, 19, 24, 29, 34, 39, 41. The dashed line is the 
theoretical values of the ACDM model. 



SNe la data dominate, which is slightly later but more 
stringent than that reconstructed from the ESSENCE 
supernova data (z ~ 0.55 — 0.73) (55[. However, there 
is another transition redshift at which the expansion of 
the universe from acceleration to deceleration, zt — 3, 
at high redshift where GRB data dominate. However, 
this transition redshift has large 1-a uncertainties (2.4 < 
zt < 3.5) compared to the former one 0.35 < zt < 0.41, 
similar to that reconstructed from the GRB data using 
the Amati relation [49]. The reconstructed EoS of dark 
energy w{z) with the 1-a uncertainties are shown in Fig. 
9. We can find that there is a singular point at z < 5 
from the reconstructed w(z). This phenomenon has been 
known to happen and just displays the fact that w(z) 
[Eq 34] is nothing more than an effective parameter that 
(in this case) fails to describe the system correctly. On 
the other hand, the singular point may be caused by the 
absence of data in the Hubble diagram at 5 < z < 6, i.e., 
between the GRB data at z < 5 and the only two GRBs 
at z > 6 as shown in Fig. 1. 

From Fig. 7-9, we can read that Hq = 
66 kms _1 Mpc _1 , q = -0.82 and w = -1.19. We can 
also find some features in the reconstructed H(z), q(z), 
and w(z), which seem to favor the oscillating models |85l - 
However, it is noted that because the available data 
of GRBs at high redshift are still quite rare now, and 
their statistical scatters are relatively large compared to 
that of SNe la, these tentative results might be artifacts, 
but nevertheless deserve further investigations. 
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FIG. 9: The reconstructed EoS of dark energy w(z) with the 
1-cr uncertainties. The red line indicates the best-fit value 
of the reconstructed w(z) which is obtained at n = 16 and 
the blue lines represent 1-<t uncertainties at n = 9,41. The 
dashed line is the theoretical values of the ACDM model. 



V. SUMMARY AND DISCUSSION 

Because of the lack of enough low red-shift GRBs to 
calibrate the luminosity relation, GRBs could not be used 
reliably and extensively in cosmology until recently. In 
our previous paper, we have proposed a new method 
to calibrate GRB luminosity relations in a completely 
cosmology-independent manner to avoid the well-known 
circularity problem [45| . 

In this work, with the recent GRB data at high rcdshift 
whose distance moduli are calibrated with the interpolat- 
ing method [45[ from the Union set of 307 SNe la folj . 
as well as the Constitution set of SNe la (6f|,the CMB 
observation from the WMAP5 result (R, l a , n b h 2 ) [6r| . 
the BAO observation from the spectroscopic SDSS DR7 



the x-ray baryon mass 
69] , and the observational 



galaxy sample (do. 2, do. 35) |£8 
fraction in clusters of galaxies ^_ 
H(z) data [z3[ll|, we find that the ACDM model is con- 
sistent with the joint data in l-o~ confidence region; this 
confirms the conclusion of many previous investigations. 
We also find that the current GRB data are substantially 
less accurate than the SNe la data, if each data set is used 
alone, consistent with previous investigations. 

The new results and insights we have obtained in this 
work are briefly summarized as follows: 

(1) In order to combine GRB data into the joint obser- 
vational data analysis to constrain cosmological models, 
we provide a simple method to avoid any correlation be- 
tween the SNe la data and the GRB data; in this method 
those SNe la data points used for calibrating the GRB 
data are not used. 



(2) Comparing to the joint constraints with GRBs and 
without GRBs, we find that the contribution of GRBs to 
the joint cosmological constraints is a slight shift in the 
confidence regions of cosmological parameters to better 
enclose the ACDM model. 

(3) Finally we reconstruct the Hubble parameter H(z), 
the deceleration parameter q(z), and the EoS of dark en- 
ergy w(z) of the acceleration Universe up to z > 6 with 
the distance moduli of the Constitution set of SNe and 
GRBs and find some features that seem to favor oscil- 
latory cosmology models; however further investigations 
are needed to better understand the situation. 

For considering the use of GRBs for cosmology, the 
gravitational lensing effect may need to be considered 
[88j ]. However Schaefer found that the gravitational bi- 
ases of GRBs are small fl7J. Recently, some possible 
observational selection bias |89h92| | and evolution effects 
[931 - 4961 ] in GRB relations have also been discussed. How- 
ever, Ghirlanda et al. confirmed the spectral-energy re- 
lations of GRBs observed by Swift [93]. Moreover, it is 
found that instrumental selection effects do not dominate 
for the Amati relation [98| and for the -E pG ak — L relation 
[99l |. as well as no strong evoluti on w ith redshift of the 
Amati relation can be found (98l. llOdt ] - Nevertheless, for 
considering GRBs as standard candles to constrain cos- 
mology, further examinations of possible selection bias 
and evolution effects should be required in a larger GRB 
sample. 

Different dark energy models may have very different 
Hubble diagrams at high redshifts [TtJ ■ Therefore, the 
best plan to investigate the property of dark energy is 
measuring the dark energy over a wide range of red- 
shifts. GRBs can extend the Hubble diagram to much 
higher redshifts beyond SNe la data [l7|. It is worth 
noticing that GRBs are important potential probes for 
cosmic history up to z > 6. However, the contribution 
of GRBs to the cosmological constraints would not be 
sufficiently significant at present, due to the large sta- 
tistical scatters of the relations and the small dataset of 
GRBs. As discussed in [ljj, GRBs are almost immune 
to dust extinction, whereas in the case of SNe la obser- 
vations, there is extinction from the interstellar medium 
when optical photons propagate towards us. On the other 
hand, SNe la are substantially more accurate standard 
candles than GRBs, which will lead to tight constraints 
on cosmological parameters; whereas a single GRB at 
high redshift will provide more information than a single 
maximal redshift SNe la 17]. Through Monte Carlo sim- 
ulations, future prospect of probing dark energy parame- 
ters with a larger sample of GRBs has been investigated 
[M HIL HH HII • It has been found that cosmological con- 
straints would improve substantially with more simulated 
GRBs expected by future observations. Recently, Xiao 
and Schaefer have compiled 107 long GRBs with their 
spectroscopic/photometric redshifts measured |l0l| . ob- 
served by BATSE, Konus, HETE, and Swift. Along with 
more and more GRBs observed from Fermi Gamma-ray 
Space Telescope (formerly known as GLAST) with much 
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smaller scatters, and its combination with the increasing 
Swift data, GRBs could be used as an additional choice 
to set tighter constraints on cosmological parameters of 
dark energy models. 

Appendix A: The SN la data used in the 
interpolating procedure for calibrating the GRB 
data 

Here we list the 40 SNe la data points that have been 
used to obtain the distance moduli of the 27 GRBs at z < 
1.4 in the interpolating procedure, and have consequently 
been excluded from the SNe la sample used to the joint 
constraints. 
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